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We propose an efficient method to filter out single atoms from trapped ensembles with unknown 
number of atoms. The method employs stimulated adiabatic passage to reversibly transfer a single 
atom to the Rydberg state which blocks subsequent Rydberg excitation of all the other atoms within 
the ensemble. This triggers the excitation of Rydberg blockaded atoms to short lived intermediate 
states and their subsequent decay to untrapped states. Using an auxiliary microwave field to care¬ 
fully engineer the dissipation, we obtain a nearly deterministic single-atom source. Our method is 
applicable to small atomic ensembles in individual microtraps and in lattice arrays. 

PACS numbers: 32.80.Ee, 32.80.Qk, 37.10.Jk, 37.10.Gh 


Introduction .— Deterministic preparation of single 
trapped atoms, or arrays thereof, is an important task 
required for many potentially interesting applications, in¬ 
cluding quantum computation and simulation Op¬ 

tical lattice experiments can realize Mott insulators with 
single atom per lattice site [§] with low particle or hole 
defect probability @0. Achieving single-atom occupa¬ 
tion of spatially separated microtraps with high fidelity is 
more challenging. Several techniques have been explored 
to reduce the trap occupation number fluctuations m- 
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Here we propose a method to filter out single atoms 
from trapped ensembles with unknown number of atoms. 
Our scheme relies on the strong, long-range interaction 
between atoms in high-lying Rydberg states [HI, HH 
which blocks resonant optical excitation of more than 
one Rydberg atom within a certain interatomic distance 
S®. We use stimulated adiabatic passage tech¬ 
nique (STIRAP) [34} to transfer any one atom of the 
ensemble from the ground state to the Rydberg state 
by a two-photon transition via an intermediate excited 
state |35H38| . This atom then strongly suppresses Ryd¬ 
berg excitation of the other atoms due to the interaction- 
induced level shifts. The Rydberg-blockaded atoms now 
turn into strongly-driven two-level atoms which rapidly 
decay from the (intermediate) excited state back to the 
initial trapped ground state [36| as well as to other un¬ 
trapped states. After several excitation and decay cycles, 
all of the Rydberg-blockaded atoms are removed from 
the trap. The reverse adiabatic transfer of the remaining 
Rydberg atom back to the initial ground state completes 
the preparation of a single trapped atom. We achieve 
high efficiency and reliability of the single-atom source 
by carefully engineering the atomic excitation and dissi¬ 
pation. The former is controlled by the time-dependent 
amplitude and detuning of the laser field, while the lat¬ 
ter can be accomplished by employing an auxiliary mi¬ 
crowave field to open appropriate decay channels to the 
untrapped states. 


Stimulated adiabatic passage with a single-atom. — For 
what follows, it will be helpful to recall the essence of 
adiabatic transfer (STIRAP) of population in an isolated 
three-level atom using a pair of coherent laser pulses [Hj] . 
A laser field with Rabi frequency UL ge couples the stable 
ground state | g) to an unstable (decaying) excited state 
|e), which in turn is coupled to another stable state |r) by 
the second laser field of Rabi frequency f l er [Fig. Oja)]. 
In the frame rotating with the frequencies of the laser 
fields, the atom-field interaction Hamiltonian reads 

V af = hA e |e)(e| + h(n ge \e)(g\ +D er |r)(e| +H.c.), 

where A e is the detuning of the intermediate state |e) 
and we assume two-photon resonance between \g) and 
|r). The eigenvalues of V a f /Ti are Ao = 0 and A± = 

A e /2 ± y/fl ge + + (A e /2) 2 and the corresponding 

eigenstates are given by 

IV'o) = (^er | g) - |r})/w 0 , 

|V’±) = (ttge I g) + A± |e) + Uler \r))/U±, 

where Wq = fl ge + Ul 2 r and uj± = + \\. The zero- 

energy “dark” state |-0o) does not contain the rapidly 
decaying state |e), which, however, contributes to the 
energy-shifted “bright” states \if±) making them unsta¬ 
ble against spontaneous decay. By adiabatically chang¬ 
ing the dark state superposition, the STIRAP process 
can completely transfer the population between the two 
stable states |g) and |r) without populating the unstable 
state |e). To this end, with the atom initially in state 
| g), one first applies the Q er field (D er . D ge ), resulting 
in (g\ i/jo) = 1. This is then followed by switching on the 
Ul ge field, leading to (r|'0o) = — D ge /wo, which results in 
|(r|?/>o)| — 1 when Q ge Q er . For sufficiently smooth 
switching, dtuio -C wo|^± — Ao|, the atom adiabatically 
follows the dark state \ipo), and the bright states 
and thereby |e), are never populated. The reverse pro¬ 
cess, at the end of which fl er fl ge , can transfer the 
atom in state |r) back to state \g). 
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FIG. 1. Level scheme of atoms driven by a pair of laser 
fields with Rabi frequencies H ge and f l er on the transitions 
|g) —» \e) and |e) —> |r). The atoms in the ground state |g) 
are confined within a volume of linear dimension L ~ 2 /im 
smaller than the blockade distance db determined by the in¬ 
teraction Vaa between their Rydberg states |r) = nSi /2 with 
large principal quantum number n (~ 50) and small decay F re 
(~ 3 KHz) and dephasing 7 r (< 0.1 MHz) rates. In (a) the 
ground state of 87 Rb atoms is | g) = 5 <Si /2 |F = l,m,F = 1) 
while the intermediate state |e) = 5 Pi /2 | F = 2,rriF = 2} de¬ 
cays back to | g) with the rate T eg = |r e (F e = 36 MHz) and 
to other untrapped states |s) = 55 i /2 \F = 2, rriF = 1, 2) with 
the rate T es = Ir e (branching ratio r es /r eg = 1). In (b) the 
transition between the ground | g) = 5 Si /2 |F = 2,rriF = 2) 
and intermediate |e) = 5 P 3/2 |F = 3, rar = 3) states is closed 
(|e) decays only to | g) with rate T e = 38 MHz), but a mi¬ 
crowave field coupling |e) to, e.g., |e') = 5 Pi / 2 |F = 2, raj? = 
2) with Rabi frequency fl ee , can open a weak decay channel 

r es = 8| ^;' 2 ~ 2.4 MHz (« F eg = T e + 4|n r ee ; 12 ) to the un¬ 
trapped states |s) = 5Si/ 2 \F = 1, 2, mp = 1). The hyperfine 
and/or Zeeman shifts of states |s) suppress their coupling to 
the H ge field, while facilitating selective interaction with a 
pushing laser of appropriate frequency and polarization. 


The multiatom system. — Consider now an ensemble of 
N atoms initially in the ground state | g) confined within 
a small trap. All the atoms uniformly interact with two 
laser fields on the transitions | g) —I |e) and |e) — 1 |r) 
with Rabi frequencies Q ge and f2 er , as shown in Fig.[lja). 
For each atom j, the interaction Hamiltonian reads 

V af = hA eV J ee + K^g^lg + M er a 3 re + H.C.), (1) 

where a 3 ^ = \p)jj{v\ are the atomic operators. The 
intermediate excited state |e), possibly detuned by A e , 
decays with the rate T eg to the trapped ground state \g), 
and with the rate r es to another untrapped lower state 
|a) which is decoupled from the driving lasers. For com¬ 
pleteness, we also take into account the typically much 
smaller decay rate T re of the highly excited Rydberg state 
|r) and its possible dephasing y r (with respect to the 
other states) due to nonradiative collisions, inhomoge¬ 
neous trapping potential, external field noise, laser phase 
fluctuations or Doppler shifts. The corresponding Lind- 
blad generators for the decay and dephasing processes 
are L 3 ge = y/F ~d 3 ge , L 3 se = y/f7sHei Hr = VHld J er , and 
H = V^v ( 2 ° 3 rr - 1 J ), where V = cr 3 gg + a 3 ss + a 3 ee + a 3 rr . 


The long-range atom-atom interactions are described 
by the Hamiltonian 

Ki = h A (dij)a l rr . <g> a 3 rr , (2) 

where A (dij) = C e /d^ g is the van der Waals potential 
between pairs of atoms i, j in the Rydberg state |r) sepa¬ 
rated by distance dij. When the interaction-induced level 
shifts A (dij) of states 1 7777 ) exceed the two-photon ex¬ 
citation linewidth of the Rydberg state w ~ 

(r e = F eg +r es ), double (multiple) Rydberg excitations 
are blocked. We can then define the blockade distance via 
A(db) X) Wmax leading to d h = ^/C 6 /w; max . To achieve 
nearly complete blockade of more than one Rydberg exci¬ 
tation, we assume that all the atoms are initially trapped 
within a volume of linear dimension L < |db correspond¬ 
ing to A (dij) > 10w max V i,j e [1 ,N]. 

Numerical simulations. — We simulate the dissipative 
dynamics of the system of N atoms usin g th e quantum 
stochastic (Monte Carlo) wavefunctions [39j, which al¬ 
lows us to treat nontrivial numbers N < 10 of strongly- 
interacting atoms exactly, without truncating the cor¬ 
respondingly large Hilbert space. In such simulation, 
the state of the system Ilk) evolves according to the 
Schrodinger equation d t |T) = — |T) with an effec¬ 

tive Hamiltonian 

n = u- l -hL 2 , (3) 


n-L+n-b 

_ ge 1 er 

^ 20 ^+ 11/4 


where 


^ = E V af + E V - 

j i<0 

is the usual (Hermitian) Hamiltonian, while 



j 


= ^[(r eg + v es )ai e + r re a 3 rr + lT V] 

j 

is the non-Hermitian part which does not preserve the 
norm of |T). The evolution is interrupted by random 
quantum jumps IT) — > L 3 gv |T) of individual atoms with 
probabilities determined by the corresponding weights 
(T| L 3 } v L 3 lv |T). In a single quantum trajectory, the 
normalized wavefunction of the system at any time t 
is given by |T(f)) = ^it))/^ (T(t)|T(t)). The ex¬ 
pectation value of any observable O of the system is 
then obtained by averaging over many, M y>> 1, in¬ 
dependently simulated trajectories, ( O) = Tr [p<D] = 
JJ Sm (^rnl O |T m ), while the density operator is given 

by Pit) = A YZ, l^m(*))(^m(0l • 

Starting from all the atoms in the trapped ground 
state, l’F(O)) = |<?i, <72 , - - -, <?jv) , we assume that during 
the evolution the atoms that decay to state |s) decou¬ 
ple from the laser field and are eventually lost from the 
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trap (which can be accelerated by a pushing laser). The 
probabilities Pjv(n) = <n£ } ) for n, out of the initial TV, 
atoms surviving in the trap are then given by the projec¬ 
tors n® EE n£i Ks, n$ = Ef = 1 (i j - 
etc. The mean number of trapped atoms is then (n) = 

E„nPjv(n) = (E^liClP - a j aa )). Typically, the initial 
number of trapped atoms is uncertain, with the Poisson 
distribution 


Pinit(N) = 


R N e -N 

N\ 


around some mean N. Then, the probabilities P(n) 
for n atoms surviving in the trap are given by P(n) = 
EN P init(N)P N (n). 

We have performed numerical simulations for various 
initial number N of atoms in the trap, results of which 
are shown in Figs. [2][4] In all simulations, we first ap¬ 
ply the field coupling states |e) and | r) with the Rabi 
frequency Q er which stays constant throughout the pro¬ 
cess, and then send a smooth pulse driving transition 
|< 7 ) —> |e) with the Rabi frequency fl ge (f) whose peak 
value is much larger than f l er . We thus expect that a 
single atom in the trap will be adiabatically transferred 
from the ground state |g) to the Rydberg state |r) and 
block Rydberg excitations of all the other atoms initiat¬ 
ing thereby their optical pumping through the interme¬ 
diate excited state |e) to the decoupled and untrapped 
state(s) |s). At the end of the process, when fl ge is re¬ 
duced to zero, the remaining Rydberg-state atom will be 
transferred back to the trapped ground state | g). 

Results. — Consider first the resonant excitation, A e = 
0, of the atoms with the fl ge and Q er fields shown in 
the upper panel of Fig. [2] The branching ratio for the 
decay of the atoms from the intermediate excited state 
\e) to the untrapped state |s) and back to |g) is large 
r es /r eg = 1, as per Fig. Ufa), while we neglect for the 
moment the dephasing of the Rydberg state, = 0. We 
find that even before the Q ge pulse attains its maximum, 
most of the atoms are optically pumped into the un¬ 
trapped state(s) |s), while the mean number of atoms in 
the trap drops below unity: P(0) ~ 0.3 and P(l) ~ 0.7, 
while P(n > 1) ~ 0. The reason for this unexpected 
and unfortunate behavior is that each quantum jump of 
a blocked atom from state \e) to either l<A or |s) projects 
a single atoms to the Rydberg state |r) [3(J whose overlap 
with the dark state |(^ok)| 2 = ^ ge /(^ge + ^er) ' s smaller 
than unity, unless Q ge fl er . Since for Q ge < f l er 
the Rydberg atom is not in the dark eigenstate, but 
has large overlap with the bright eigenstates as well, 
|(^±|r)| 2 = fl 2 r /(2Q ge + 2Ql r ), it is coupled to the state 
|e) from where it can decay to | g) or |s). In a closed 
three-level system [36J, the decay to |g) would simply 
restart the excitation process of this or any other atom 
by the fields Q ge and Q er . In the present case, however, 
the decay of all the atoms but one to the decoupled state 
|s) may leave the last atom in the Rydberg state |r) be¬ 
fore the f2 ge (i) pulse reaches its peak value Q ge f l er . 




FIG. 2. Dynamics of the system, as obtained from M = 
200 — 500 independent realizations (trajectories) of the nu¬ 
merical experiment for each number N of atoms with the 
level scheme of Fig. 01 a), irradiated with two resonant fields, 
A e = 0, of Rabi frequencies Q, ge and Q, er shown in the upper 
panel. Main panel shows the mean number (n) of surviving 
atoms in the trap, starting with N = 1 — 10 (from bottom to 
top) ground state atoms, for ■y r = 0. Inset shows the corre¬ 
sponding probabilities Pn{u) of finding n trapped atoms in 
state | g) at the end of the process (open bars), while the black 
('Yr = 0) and gray ( 7 ,. ^ 0) circles denote the P(n) averaged 
over the initial Poisson number distribution Punt (TV) of the 
trapped atoms with the mean N = 5. 


The probability of loosing this last remaining atom can 
thus be estimated as 


r es + r 


-(l-|(^o|r)| 2 )^ 


r fr 

1 es* L er 


eg 


(r e 


eg 


mie 


Mir V 


which, with r es = r eg and f l ge — Q er is about 1/4, 
consistent with (n) and P(0) in Fig. [ 2 ] 

Including now in the simulations realistic dephasing, 
7 r 0 , amounts to additional coherence relaxation be¬ 
tween the long-lived states \g) and |r) with the rate 
7 rg = Irre+ 27 ,.. As is well known 341, the STIRAP in a 
three-level atom is very sensitive to such dephasing, since 
the decoherence of the dark-state superposition leads to 
additional population of the bright states (36], which in 
the present system can decay to the untrapped state(s) 
|s). Hence, in the inset of Fig. [2] we observe further de¬ 
crease of the survival probability P(l) ~ 0.63 for a single 
trapped atom. 

To increase the single atom survival probability P(l), 
we should therefore postpone the decay of atoms from 
state |e) to |s) until f l ge is much larger than f l er . We 
thus take a finite detuning A e of intermediate level |e) 
and switch-on S2 ge slightly faster. The results of simula¬ 
tions are shown in Fig. [3] were we indeed notice signif¬ 
icant improvement of the single-atom filtering, P(l) ~ 
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FIG. 3. Same as in Fig. [21 but with finite detuning A e and 
slightly different optimized pulse shape of ft ge . 


0.87 or 0.8 with additional dephasing ^ r . These results 
are obtained after some optimization of both f2 ge and A e , 
but for the given atomic level scheme further substantial 
improvements do not seem possible. 

To achieve nearly perfect filtering, we now consider the 
closed three-level system of Fig.OJb), where state |e) de¬ 
cays spontaneously only to the trapped state |g), while 
we can engineer an additional decay channel to the un¬ 
trapped states |s). To this end, we use a microwave field 
which weakly couples |e) to an auxiliary state le') with 
the Rabi frequency Q ee > that is small compared to the 
total decay rate T e / of |e / ). Upon adiabatic elimination 
of |e'), we obtain a small additional decay on the tran¬ 
sition |e) —► |g) with rate 4|U ee /| 2 /3T e /, and a decay 
to the untrapped states with rate T es = 81 O ee ' ] 2 /3r e ' 
which can be controlled via the microwave field inten¬ 
sity |f 2 ee /| 2 (the numerical prefactors are related to the 
branching ratio of the decay of state |e') to | g) and |s), 
r e 'g/r e ' s = 1/2). FigureS] shows the results of our simu¬ 
lations demonstrating dramatically improved single atom 
filtering, P(l) ~ 0.96 [P(0) ~ 0.02 and P( 2 ) ~ 0.01, 
P(3) ~2x 10 -4 ], at the expense of a longer duration of 
the process. Even sizable dephasing 7 r of the Rydberg 
state with respect to the ground state does not signifi¬ 
cantly reduce the efficiency of the single-atom filtering, 
P(l) ~ 0.92 [P(0) ~ 0.05], which is due to the robustness 
of adiabatic passage in a Rydberg-blockaded ensemble 
composed of atoms with nearly-closed transitions [36j ■ 

We finally note that in obtaining the averaged proba¬ 
bilities P(n) of the final number of atoms, we assumed 
the initial Poisson distribution Pi n n(N) of the atom 
number in the trap with the mean N = 5, for which 
-Pinit(O) ^ 0.007 and X)jv>io ~ 0 . 02 , and per¬ 

formed simulations for N < 10. The initial 2% popula¬ 
tion of the trap not included in our simulation can only 


FIG. 4. Same as in Figs. [2] and 0 but for the scheme of 
Fig. mb) with resonant fields, A e = 0, and weaker decay 
to |a) (notice the longer time scale of the process). Thin 
dotted lines in the main panel show single trajectories for the 
corresponding initial number of atoms N = 1 — 10. 

increase the final single atom probability. Our method 
should work even better for larger initial number of atoms 
in the trap, which is, however, prohibitively difficult to 
simulate numerically. 

Summary .— To conclude, we have proposed and an¬ 
alyzed a viable scheme to filter out single atoms from 
mesoscopic ensembles containing unknown number of ini¬ 
tially trapped atoms. We have shown that the combi¬ 
nation of two-photon adiabatic transfer of the atoms to 
the strongly interacting Rydberg state and carefully en¬ 
gineered decay of atoms from the intermediate excited 
state can lead to nearly-deterministic production of single 
trapped atoms in individual microtraps or optical lattice 
potentials, which is an important prerequisite for quan¬ 
tum computation with atoms serving as qubits or quan¬ 
tum simulations of many-body spin systems. Moreover, 
our studies represent an important and interesting new 
facet of the very active field of dissipative generation of 
states and processes in open many-body systems [13, EH . 
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